Density Singularities and Cosmic Structures 

G. Murante, A. Provenzale, E. A. Spiegel^ and R. Thieberger* 
Istituto di Cosmogeofisica, Corso Fiume 4, 1-10133 Turin, Italy 
t normally at: Dept. of Astronomy, Columbia University, New York, NY 10027, USA 
* normally at: Dept. of Physics, Ben Gurion University, Beer Sheva, Israel 

M.N.R.A.S., accepted 

Abstract 

We show that an arrangement of singular density concentrations accounts for the scaling exponents 
observed in the luminous matter distribution in the Universe for scales below lOMpc. This model 
provides a good description of the matter distribution on those smaller scales. 

1. Introduction 

Millennia ago, people looked at the stars in the sky and saw the constellations. Today we look at 
catalogues of galaxy positions (Shane & Virtanen 1967, Zwicky et al. 1968, Maddox et al. 1990) 
and redshifts (Huchra et al. 1990, Haynes & Giovanelli 1988, da Costa et al. 1989, Shectman et al. 
1996) and see groups, clusters and superclusters (de Vaucouleurs 1970, Bachall 1988) and, on larger 
scales, voids, filaments and walls (see Jones 1992 for an introductory review). The big difference 
is in our present ability to quantify the distribution and thus give meaning to the identification of 
these structures. Nevertheless, we still retain the desire to rationalize what we see and there are 
several ways to go about this, as described by Peebles (1980, 1992) and Borgani (1995) as well as 
authors in Bonometto, Primack and Provenzale (1997). The same quantitative approaches may 
be applied to the results of simulations of gravitationally interacting point masses in an expanding 
universe, since these reveal a similar variety of structures. Agreement between the analysis of 
observations and the computer results encourages the belief that there may be an underlying, 
simple theoretical picture of the evolution of the galaxy distribution. 

We may also look for a rationalizing statistical picture and there are two that now seem to be 
the main contenders for acceptance. The first is the conventional fractal picture (Mandelbrot 1982, 



Martinez 1991, Provenzale 1991, Coleman & Pietronero 1992, Borgani 1995) whose antecedents are 
traceable to the past century. This fits in well with the belief that cosmic structures are formed in a 
cascade process like that of turbulence theory (Thieberger, Spiegel & Smith 1990, Provenzale 1991, 
Provenzale et al. 1992) with a resulting rich geometrical texture. The second, equally compelling 
view, is that local density singularities form and organize the matter distribution around them 
(Peebles 1974, McClelland & Silk 1977, Gurevich and Zybin 1988; Sheth k Jain 1996, White 
1996). 

Our aim in this paper is to scrutinize the singularity picture by testing its ability to meet the 
quantitative tests and produce scalings like those seen in the galaxy distribution. That is, we shall 
find the generalized dimensions of a set of singularities, just as one does in analyzing any fractals. 
As we shall see, the singularity picture is capable of describing the scaling properties of the galaxy 
distribution well on small scales of up to 10 Mpc. 

2. Statistical descriptors. 

Analyses of the spatial arrangement of galaxies have been based mainly on both the correlation 
function (e.g., Peebles 1980) and the correlation integral (Grassberger & Procaccia 1983, Paladin 
& Vulpiani 1984). The two-point correlation function, £, of a given set of galaxies is defined such 
that the joint probability of finding a galaxy in a volume element 5V\ and another in a volume 
element 5V2 a distance r away from the first is given by the expression 5 2 P = n 2 SViSV2[l + £(r)], 
where n is the mean number density of the sample. 

For a uniform Poisson distribution, £ = and the galaxies are uncorrelated. If £ is positive, 
the galaxies are said to be correlated. In the latter case, as r goes to zero, £ becomes large like an 
inverse power of r. This behaviour is expressed in the form £ oc r~ 7 where 7 is a constant. The 
precise value of 7 depends on the morphological type of galaxies that are considered and possibly 
on their luminosities. The value 7 « 1.8 has been found for r less than five or ten h~ x Mpc, 
especially from angular position catalogues (Peebles 1980), where h is the Hubble constant in units 
of 100 km/s/Mpc. 

Higher order correlations may also be calculated, though the formalism for such richer de- 
scriptions has been more actively developed for the generalised correlation integrals. In recent 
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elaborations of the theory of point sets (see the reviews by Paladin & Vulpiani 1987 and by Bor- 
gani 1995), one introduces the moments 



where q is a parameter and is the Heaviside function. The inner summation is over the whole 
set of N — 1 galaxies with coordinates Xj, j ^ i, and the outer summation is over a subset of N' 
galaxies, taken as centers, with coordinates Xj. 

To get an intuitive feel for this formula, consider M(r), the average number of galaxies of the 
set lying within a distance r of a typical galaxy. We see from the formula that C^r) = J\f(r)/N. 
We may approximate J\f(r) as the integral, 4tt J Q r n(s)s 2 ds, where n(r) is the density function. 
The existence of structure in the distribution of galaxies is connected to the dependence of n on 
r and it may be shown that for a suitable constant nh, we may write n(r)/nh = 1 + £(r) where 
£(r) is the pair correlation function just defined (Provenzale 1991). In the last few years, many 
lively discussions have focussed on the question whether rih becomes independent of the sample 
size above a certain threshold (see e.g. Lemson & Sanders 1991, Coleman & Pietronero 1992, 
Provenzale, Guzzo & Murante 1994). 

More generally, for any integer value of q, C q (r) is the fraction of (/-tuples in the set whose 
members lie within a distance r of one another. For sufficienly small r, C q will go to zero and, 
for a typically well-behaved set, it will vanish like r Dq , where the index D q is called a generalized 
or Renyi dimension (Renyi 1970, Halsey et al. 1986). The quantity C2 is called the correlation 
integral, since it is the integral of 1 + £. Consequently, by looking at small r, we learn that 
D 2 = 3 — 7 for those scales where 1 + £(r) may be approximated by £(r). 

Since whatever scaling exists cannot be expected to be perfect, the values obtained for D q 
(and 7) depend on the range of scales of r under consideration. To simplify the discussion, one 
may at the outset limit the range of r to include only portions of the range with a single slope. 
However, when the existence of a scaling range is not guaranteed a priori, the practice of using 
a local slope is often followed. In this method, one plots logC q against logr and determines the 
local slope D q (r) as a function of the scale r. In any range of r for which D q (r) is nearly constant 
there is said to be scaling behaviour for that particular range. 




(1) 
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A certain amount of ambiguity in interpreting the observed galaxy distribution may arise 
because some investigators try fit a power law to £ while others prefer fitting one to C2, thus to 
1 + £. Evidently, for small enough r, this makes no difference, but discrepant results have been 
reported because practical considerations sometimes force the range of r included in the analysis 
to be too large (Calzetti & Giavalisco 1991). Such discrepancies are especially noticable in denser 
sets with strong self-overlap (Thieberger et al. 1990), even in theoretical sets whose properties 
are well controlled. For real data, failure to restrict to very small values of r may aggravate the 
disagreements. Nevertheless the earliest results giving 7 fz 1.8 seem robust for small enough r. It 
is a challenge to the theory to explain this number while rationalizing the existence of the different 
scaling regimes. 

Other results of the analysis of the galaxy distribution (e.g., Martinez & Jones 1990, Borgani 
1995) and of the output from numerical simulations (Valdarnini, Borgani & Provenzale 1992, 
Colombi, Bouchet & Scheffer 1992, Yepes, Dominguez-Tenreiro & Couchman 1992, Murante et 
al. 1996), confirm the existence of approximate scaling behaviour with 1 < D q < 1.4 for q > 2, 
at scales below 5 h~ x Mpc. At larger scales, the existence of another scaling regime has been 
suggested (Guzzo et al. 1991, Yepes et al. 1992, Murante et al. 1996) with D 2 ~ 2. At still larger 
scales, a growth of the generalized dimensions with r is usually observed (Peebles 1989). This 
behaviour of the generalized dimensions at the largest values of r, with D q (r) increasing with r, is 
in general considered to reflect an approach to homogeneity at large scales. This is in keeping with 
the smoothness of the background radiation, and most workers expect the dimension to approach 
three, as for a space-filling distribution. In any case, at the largest scales, other factors like time 
dependence and curvature will have to be included before serious conclusions are reached. For the 
richer structures seen at smaller scales, these problems do not intrude and that is the regime we 
study here. 

3. The Singularity Picture. 

Numerical simulations of structure formation are carried out in an expanding background typically 
beginning at a state that corresponds to the linearly evolved density field produced at the time 
of recombination. The initial perturbations to the uniform state are often made as large as the 
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observed fluctuations in the background radiation allow. The subsequent evolution shows the 
formation of mass concentrations of various kinds — flat structures, filaments and clusters - 
and generically leads to the development of a number of local strong condensations (Efstahiou et 
al. 1988, Zurek and Warren 1994). These structures form in the post-recombination era from 
density perturbations that have survived either viscous and radiative damping and dispersion due 
to random velocities. 

The largest of the evolving structures are relatively immune to the effects of pressure after 
recombination. Any deviations from sphericity in them will be exaggerated (see the literature cited 
in Peebles 1980) as these early formations take on the flattened form known as Zeldovich (1970) 
pancakes. These arrange themselves in a complex skein of sheet-like structures according to many 
simulations of this process. At the same time, as long as the pressure remains unimportant, the 
collapse will continue to form secondary structures, particularly at the intersections of pancakes. 
As many discussions of the limiting behaviour of these structures suggest, the density will locally 
evolve toward the formation of singularities (Efstahiou et al. 1988; Gurevich and Zybin 1988; 
Navarro, Frenk & White 1995; Cole & Lacey 1996; White 1996). 

The creation of local clusters seen in the real and simulated galaxy distributions has focussed 
attention on spherical collapse. This leads to the formation of singular densities in finite times 
with a density n(R) oc R~ a where R is the distance from the singularity. Singularities of this form, 
with a between 1.5 and 3, are self-similar solutions of the Vlasov-Poisson equations (Henriksen & 
Widrow 1995). Calculations by Gurevich & Zybin (1988) for special, but plausible, pressure-free 
initial conditions lead to a = 24/13. The case a = 2 corresponds to an isothermal gas sphere, with 
pressure. Tight spherical structures do seem to be the outcome of many simulations of structure 
formation (Peebles 1980, White 1996) though there is no clear argument for the existence of such 
spherical, presumably nongeneric, structures. 

We mention these theoretical notions as motivational and do not assess their validity here. Nor 
do we intend to compare their importance to those of the other processes that have been considered 
in discussions of structure formation, such as cosmic strings (e.g., Kolb & Turner 1990) and fractal 
cascades (Provenzale, 1991). For while it is true that the power-law behaviour of the correlation 
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integrals and the fractional values of the generalized dimensions may reasonably be attributed 
to a fractal nature of the galaxy distribution (Efstathiou, Fall & Hogan 1979, Mandelbrot 1982, 
Coleman, Pietronero & Sanders 1988, Provenzale 1991, Martinez 1991, Borgani 1995), we are 
interested here in the singularity picture as a fruitful way to think about the galaxy statistics. 

4. The Basic Scalings. 

We now show that the scaling behaviour of a random superposition of singularities of a given type 
provides a good description of the higher order statistical moments found in the galaxy distribution. 
These results extend previous calculations of the shape of the correlation function £ for an isolated 
singularity (Peebles 1974) or for a distribution of density singularities (McClelland & Silk 1977, 
Sheth & Jain 1996). 

Consider first the generalized dimensions of the set of points disposed around a single con- 
densation with power-law profile. For large, positive q, the correlation integrals are dominated by 
the very dense regions of the distribution near to the singularity itself. In the formula for C q , the 
sum over j gives the number of galaxies within a distance r of the ith galaxy. This quantity will 
clearly peak sharply for a galaxy at (or very near to) the singularity itself where it becomes equal 
to the mass interior to a sphere of radius r around the singularity. The second sum averages this 
quantity over the galaxies and it is completely dominated by galaxies at the singularity. Hence, 
for small r, with r rj R, the expression for C q (r)/[r 3 n(r)] approaches a constant and does so most 
rapidly for q — > oo. (Here R is the distance from the singularity and r is the separation of points 
in the statistical moment). 

This argument shows that a singular density distribution, with n(R) cx R~ a for small R, has a 
scaling behaviour (at least for small r) with = 3 — a, in analogy with the more familiar formula 
D 2 = 3 — 7. For any finite q, the value of the dimension D q is larger than 3 — a, because of the 
'dressing' of the dimension by the contribution of points away from the singularity. Nevertheless, 
for sufficiently strong singularities (a > 1.5), we find a good scaling behaviour even for moderate 
values of q. In particular, using the result 7 = 2a — 3 of Peebles (1974), one has that at small r 
the value of D 2 tends to D 2 = 3 — 7 = 6 — 2a on those scales where £ « 1 + £. In the singularity 
picture, the scaling exponents, D q , are thus largely determined by the singularity exponent a. 
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More explicitly, let points be put down in a three-dimensional space at random with a prob- 
ability proportional to the density distribution n(R) = no(R/R max )~ a with R < R max and with 
n(R) = for R > R ma x- To have finite masses in such structures, we adopt finite values of 
Rmax- In calculating the correlation integrals to find D q (r) for this finite sample of points, we 
face the difficulty that the formula for C q (r) is not well adapted to a set of galaxies with members 
near the edge of the set. So, in practise, we base the evaluations on only a subset of the points 
whose members are well away from the edges of the sample. We then treat the remaining points 
as background objects to complete the evaluations. We do this for a variety of such subsets and 
average over the results from each. Using this approach, we obtain the results shown in Fig. 1 
for the singular distributions discussed above, for the cases a = 1.8, 2 and 2.4. Scaling is clearly 
evident for the moments with q > 2. For q = 2, at small r the value of D 2 tends to D 2 = 6 — 2a, 
as expected from the theory. For higher values of q, the dimensions become smaller and tend to 
Doo = 3 — a for very large q. 

For an arrangement of several such singular density concentrations, the situation is like that 
in Fig. 2, which shows a random, uncorrelated distribution of singularities with a = 2. The 
number density of singularities is here n s = 10 and, for simplicity, we work in a volume of unit 

— 1/3 

size. The average distance between two singularities is d = n s ■ The number of points around 
each singularity is N = 100, 000 and we have chosen R max = 0.8d, so that different singularities 
only mildly overlap. We state the dimensions for this case here and return in the next discussion 
to the question of the effects of varying the parameters. 

In Fig. 3 we show the generalized dimensions calculated for the distribution of points of Fig. 2. 
Good scaling behaviour is found, similar to that found in the analysis of the real galaxy distribution 
(Martinez & Jones 1990, Borgani 1995) as well as in N-body simulations (Valdarnini et al. 1992, 
Murante et al. 1996). When the singularities are placed more deliberately, for instance on the 
sheet-like structures produced by large-scale pancake formation, the qualitative appearance of the 
distribution becomes even more plausible. But an extensive study of this latter issue is deferred to 
future work. Before turning to the meaning of these results we next examine how they are affected 
when the parameters are assigned different values. 
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5. Role of the parameters 

The form of the singularities discussed in the previous section is n(R) = rio(R/ R max )~ a . For 
a distribution of such singularities we have to specify the values of four parameters: (1) the 
number density of singularities n s , or, equivalently, the average distance between two singularities, 

— 1/3 

d = n s ; (2) the exponent in the power-law, a; (3) the cutoff radius, R m ax, around a single density 
enhancement; and (4) the density scale, n = n(R max ). We also have N = 4-7T J Rmax n(r)r 2 dr, 
the number of points around each singularity, which, in terms of the other parameters, is N = 
47rn i?^/(3-a). 

Clearly, R ma x and d are not independently meaningful; the relevant parameter is a = R max /d, 
which measures how much overlap the singularities have with each other. When the singularity 
centers are too closely packed together, the overlap becomes strong and it becomes increasingly 
hard to detect scaling behaviour unless the number of points in the sample is large enough to allow 
resolution of the cores of the individual density concentrations. This problem is common in dense 
sets (Thieberger et al. 1990) and needs to be borne in mind as we examine the role of the number 
of points N per singularity. In Figure 4 we show the results of the scaling analysis for the same 
distribution discussed in the previous section (n s = 10, a = 2, a = 0.8), but with iVo = 25,000 
(panel a) and N = 7,000 (panel b). For the different values of No, no difference is observed 
at large scales, while at small scales the lack of statistics destroys the scaling behaviour. This 
shows that varying N (or n ) for fixed a, a and n s does not influence the scale of the transition 
to homogeneity but only the detectability of the small-scale scaling regime. If N is too small, 
there is the danger that the small-scale regime is completely lost owing to the lack of statistics. 
Analogous results have been discussed by Borgani et al. (1993) in the case of fractal distributions. 

Next we turn to the role of the power-law exponent, a. Figure 5 shows the results of the scaling 
analysis for a distribution of singularities with n s = 10 and a uniformly distributed between 1.6 
and 2.4. Here N = 100,000 and a = 0.8. For low values of q, the scaling is less good due 
to the averaging over contributions from different singularities. For larger g's, the correlation 
integrals give more and more weight to the most overdense regions, that are associated with the 
stronger singularities. At these values of q, scaling is observed with a dimension close to that of 

8 



the singularity with the largest values of a. 

The measure of overlap between concentrations is a = R ma x/d and, for galaxies, this will 
be determined by how long the density perturbations have been evolving. If a is quite small, the 
singularities are isolated from each other and wide empty spaces remain between them. In this case, 
the scaling analysis indicates a lack of points at intermediate scales. On the other hand, if a is too 
large many points are in the overlapping portions with a corresponding lack of statistics at small 
scales. In this case, very large values of iV"o are needed in order to detect the small-scale regime. A 
value a = 0.8 provides very little overlap without leaving great empty spaces. However, the precise 
value of a is not crucial to the scaling results. Figure 6 shows such results for a distribution of 
density singularities with n s = 10, a = 2, N = 100,000 and R ma x uniformly distributed between 
0.5d and 1.5d. The results are very similar to those obtained with R max = 0.8d. 

Finally, we note that some investigations of the results of cosmological N-body simulations 
have indicated that large density peaks could produce behaviour other than pure power laws. In 
fact, if the time of singularity has not been achieved, there will be density maxima that are rounded 
in their inner cores. One form of density maximum that has been considered is (Navarro et al, 1995) 
n(R) = {Rq / R)( \ + R/ R s )~ u , where Rq is a normalization factor and R s marks the transition from 
the R~ x behaviour at small scales to the R~ v ~ x behaviour at large scales. Typically, v = 3 for the 
Hernquist model and v = 2 for the NFW model (Navarro et al, 1995). Our analysis shows that for 
an isolated density peak with this shape, no clear scaling behaviour is observed. Independently of 
whether we understand the origin of power law behaviour in the density, it does seem to be in line 
with the observations of large-scale structure. 

6. Conclusions 

The aim of this work was to see how well an arrangement of density singularities can account 
for the observed scaling of the galaxy distribution. We conclude that it does this as well as the 
conventional picture of a fractal distribution (in which we include the multifractal as well). Of 
course, one might argue that, since the dimensions we find are not integers, the singularity picture 
is also a fractal one and, on this matter, we would not disagree. Nevertheless, we believe that 
it is fair to claim that an explanation of the observed scaling in terms of singular densities is 
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qualitatively different from a conventional geometric fractal picture. 

On physical grounds, the singularity scenario is well justified and it is supported by N-body 
simulations that indicate the existence of strong density enhancements at small scales. On this 
view, the origin of the scaling exponents for small scales in the present data is to be found in the 
structure of the dominant singularity with the value 7 = 1.8 corresponding to a = 2.4 (Peebles, 
1974). At larger scales, the spatial distribution of singularities themselves becomes important 
(Sheth & Jain 1996). 

The existence of scaling behaviour in the cluster distribution (Bachall 1988, Borgani 1995) 
may indicate a geometrical arrangement of the small scale condensations. This could point to 
a truly fractal distribution of galaxy clusters or higher structures on large scales. On the other 
hand, this regime may result from another type of singularity, such as pancakes. Work is now in 
progress to properly address this issue in the framework of the singularity picture and to study the 
dependence of the density profile of the singularities on the initial conditions. 

But it is not the scaling distribution alone that decides which description to invoke in trying to 
learn whether large scale structure is formed mainly by cascades or by collapse onto singularities. 
There are other observational features indicating that singularities are prevalent on the smaller 
scales. Recent detection of cuspy density distributions that seem describable as isothermal spheres 
suggest singular density distributions with a « 2 (Fukushige & Makino 1996). Also nearly flat 
galactic rotation curves seem compatible with the present estimates of a. 

For a variety of disk galaxies, the rotation curves vary like r e with e of the order of ±0.1 
(Begeman et al. 1991). On the other hand, for the spherical singular density distribution we are 
considering, the rotational velocity will vary like A 2 ~ a ^ 2 . This gives us values of e = (2 — a)/2 in 
the observational range for a = 2 ± 0.2. The agreement in order of magnitude is suggestive of a 
universality of density singularities. There is much to be understood in the galaxy data, including 
the correlation of e with other properties of the galaxies, but this coincidence of values seems worth 
pondering. 

We conclude that the presence of density singularities is likely to be a common feature of the 
large scale structure of the universe and that the study of the evolution of such singularities will 
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continue to be a fruitful activity. 

We thank Dr. R.K. Sheth for his through reading of the manuscript. We acknowledge B. Dubrulle, 
J. v. Gorkom, I. Grenier and A. Refregier for pointing out relevant references. 
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Figure Captions. 

Figure 1. Effective generalized dimensions for q = 2, 4, 6 for points placed at random accordingly to 
the density n(R) = R~ a , for a = 1.8 (panel a), a = 2 (panel b) and a = 2.4 (panel c). Dimensions 
are obtained by averaging over ten resamplings of the subsets used as centers; error bars are the 
3cr deviations over the resamplings. 

Figure 2. Point distribution sampling the density field generated by a random distribution of 10 
singularities in the unit cube with size L = 1. Each singularity has a shape n(R) = n (R/ R max )~ a , 
here a = 2 and R ma x = 0.8<i where d = 10 -1 / 3 is the average distance between two singularities. 
The number of points per singularity is A^ = 100, 000. 

Figure 3. Effective generalized dimensions for q = 2, 4, 6 for the point distribution shown in Figure 
2. 

Figure 4. Effective generalized dimensions for q = 2, 4, 6 for a distribution similar to that shown 
in Figure 2, but with N Q = 25, 000 (panel a) and iV = 7, 000 (panel b). 

Figure 5. Effective generalized dimensions for q = 2, 4, 6 for a distribution of 10 singularities with 
a uniformly distributed between 1.6 and 2.4. Here R max = 0.8d and -/V = 100,000. 

Figure 6. Effective generalized dimensions for q = 2,4, 6 for a distribution of 10 singularities with 
a = 2 and R max uniformly distributed between 0.5<i and l.bd. Here N = 100,000. 
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Figure 1b 
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Figure 1c 
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Figure 3 
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Figure 4a 
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Figure 4b 
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Figure 5 
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Figure 6 
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